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A mechanism of the high temperature ferromagnetism in polymerized fullerenes is suggested. It 
is assumed that some of the Ceo molecules in the crystal become magnetically active due to spin 
and charge transfer from the paramagnetic impurities (atoms or groups), such as hydrogen, fluorine, 
hydroxyl group OH, amino group NH2, or methyl group CH3, dispersed in the fuUerene matrix. The 
exchange interaction between the spins localized on the magnetically active fullerenes is evaluated 
using ab initio calculations. The nearest neighbour and next nearest neighbour exchange interaction 
is found to be in the range 0.1 0.3 eV, that is, high enough to account for the room temperature 
ferromagnetism . 
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I. INTRODUCTION 

Recent developments in material science have resulted 
in the discovery of a novel class of magnetic organic ma- 
terials based on carbon, such as fuUerenesiiSiSi^iSiSiLSiSiiS 
and graphite. ^""^'^^ Fullerenes Cgo intercalated with the 
organic TDAE-molecule is ferromagnetic below 16 K.^'^ 
The ferromagnetic ordering below ss370 K in the 
PVDF-Ceo composite, where PVDF is polyvinylideneflu- 
oridc, (-CH2-CF2-)„, has been reported.^ More recently, 
room temperature ferromagnetism has been reported in 
fullerenes polymerized by means of photo-processes^i^ 
and by high pressure - high temperature treatmenliLSiSiifi 
and hydrofuUeriteSfi^i In the present paper, we focus our 
attention on the problem of the high temperature mag- 
netism in polymerized fullerenes. 

In accordance with the band structure calculations, 
an ideal lattice of polymerized fuUereneg^^ is not ex- 
pected to show magnetism. Some yet unidentified struc- 
tural or chemical imperfections are of crucial importance 
for the fullerenes to become magnetic. The data sug- 
gest that the type of polymerization, rhombohedralS&S 
or tetragonal,-'^'^''^^ is not of primary importance, and 
that the polymerization is necessary but not sufficient 
a condition for a high temperature ferromagnetism. For 
photo-polymerized fullerenes, the presence of oxygen is a 
prerequisite.''"'^'^ The results of the Ref. Q give the hint 
that Ceo radical adducts, Cgo^n where R originate from 
the organic polymer fragments, are responsible for the 
magnetism. Seeing that the intrinsic magnetism in poly- 
merized fullerenes remains an experimentally controver- 
sial issuei^iiSiii, it is of interest to demonstrate that the 
high-temperature magnetism is a theoretically feasible 
possibility in realistic assumptions about the materials. 

To answer the question about a microscopic origin of 
the ferromagnetism in carbon-based materials, the the- 
ory has to find the structural element of the carbon 
matrix that carries unpaired spins, and to show that 



the interaction between the spin of the units leads to 
a parallel spin alignment. So far, the theoretical devel- 
opment has been in the picture of magnetically active 
structural defects, i.e., defects with localized spins, inter- 
action of which is mediated by the magnetically passive 
carbon matrixj-'^i-'^i^^i^-'i^^i^^i^"' Defects of different com- 
plexity and topology were considered: carbon vacancies 
in graphite, ^^'^^ carbon adatoms on the graphene layer )2£ 
vacancies in the fuUerene cages—, partially opened in- 
termediate fuUerene cage structures with the zigzag-type 
edger^ the carbon tetrapods with negative Gaussian 
curvatureS^, and the special open-cage defect structure 
with the hydrogen atom bonded chemically to one of 
defect carbon atomsi24 It is yet unclear whether these 
defects are present in real fuUerene samples (see, how- 
ever, Ref. where atomic-scale defects were observed 
in graphene layers) and whether the spin-spin interaction 
is ferromagnetic and strong enough to account for the 
high temperature magnetism of polymerized fullerenes. 

As pointed out in Ref. |3j in the periodic version of 
the defective structure of Rh-Cgo proposed in the Ref.l2ll 
i.e., in Rh-Csg, the interaction between the magnetic mo- 
ments localized on the cages is not ferromagnetic. The 
defective structure model proposed in Ref. |2J does pro- 
vide a mechanism for ferromagnetic ordering, but the 
computed energy difference between ferromagnetic and 
antiferromagnetic states is too small (3 meV per cage) to 
explain the high temperature magnetism. Besides, the 
magnetically active structural defects approach is hardly 
applicable to cold-polymerized fuUerenesS*^*^ where the 
formation of Ceo cage defects is ruled out. A model, 
where magnetism arises in a system of undamaged buck- 
yballs, has been suggested in Ref. In accordance 

with the calculations, a neutral Ceo dimer turns into the 
triplet state provided the interfuUerene bond is shortened 
below the critical length 1.3 A.^^ This value is apprecia- 
bly smaller than the experimentally observed bond length 
1.58- 1.62 A, and an enormous pressure would be needed 
to decrease the intermolccular distance below the critical 
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one. It is not clear whether the model is universal in 
particular in the view of the observation that the mag- 
netic transition may take place even without any external 
pressurei2iliSi& 

In this paper we suggest a scenario where doping cre- 
ates fuUerene radical adducts CeoR, i-e., paramagnetic 
species with unpaired spins localized on fuUerenes. In 
our model, ferromagnetism occurs in a network of the 
paramagnetic species dispersed in a polymerized fuUerene 
matrix. The model is based on the computational 
observations! that the ground state of a doubly charged 
[2+2] cycloadduct dimer (Cqq)^ is triplet^ with the sin- 
glet state separated by the gap of order of 0.2 -;- 0.3 eVM: 
In other words, the spins localized on the neighboring 
C^Q ions interact ferromagnetically, and the interaction 
is strong enough to account for the high temperature fer- 
romagnetic transition in polymerized fullerenes. 

In our model we assume that some of the Cgo molecules 
in the crystal, initially spin inert and charge neutral, 
become magnetically active (paramagnetic) due to the 
charge transfer from radicals (impurities), naturally ac- 
companied by the spin transfer. Although the electron 
transfer is the actual reason for the fullerenes to acquire 
spin, we label the scenario "spin transfer" because Cgo 
acquire spin i for any direction - to or from the Cgo 
molecule - of the electron transfer. The purpose of this 
paper is a detailed study of the spin transfer model using 
first principle calculations. 

The present mechanism is best illustrated by the case 
of a paramagnetic impurity A {A = alkali metals, hy- 
drogen, fluorine, the hydroxyl group OH, amino group 
NH2, methyl group CH3 etc ) allowing for the charge- 
spin transfer reaction 

paramagnetic A + diamagnetic Cgo — >■ 
diamagnetic A'^{A^) + paramagnetic Cgo(C^Q) (1) 

with the formation of charged paramagnetic ions C^q. 
Here we speak about the final states of the charged im- 
purity or fuUerene, with particular emphasis on their spin 
states. One can see that our scenario has many features 
in common with the charge transfer complexes model'^ 
proposed for the description of magnetic properties of 
TDAE-Cgo and to the qualitative picture of paramag- 
netic (open shell) Ceo radical adducts considered respon- 
sible for the ferromagnetism in Cgo dispersed in PVDF.^ 
In the model under consideration, the primary source 
of the local spin moments is the radical impurities dis- 
persed in the fuUerene lattice. In this respect, the spin 
transfer magnetism in polymerized Cgo is akin to the fer- 
romagnetism in a dilute system of magnetic atoms in an 
insulator or semiconductor, with some important differ- 
ences however. First, the reaction Eq. transforms 
paramagnetic impurity atom (or molecule) to a diamag- 
netic ion with closed electron shells. This means that 
the impurities, nominally paramagnetic, play only a pas- 
sive role in our scenario of the sources of net charges 
and spins. Second, both magnetically active centers (ac- 
tual paramagnetic species) and magnetically passive ones 



(diamagnetic matrix) are constructed from the fuUerene 
molecules. 

To validate our physical picture of magnetism, we eval- 
uate the coupling between paramagnetic species (nearest 
neighbors and next nearest neighbors) in the polymer- 
ized fuUerene matrix. The role of the fuUerene radical 
adducts, CeoR, where i?=H, F, OH, NH2 or CH3 in the 
formation of the ferromagnetic ground state is studied. 
For this purpose, the first-principles cluster calculations 
of the electronic structure, optimized geometry, and en- 
ergies of the two low-lying levels corresponding to the 
singlet and triplet spin states for the pairs of ions CgQ 
and Ceo radical adducts, CqoR (-R=H, F, OH, NH2, CH3) 
connected by [2-1-2] cycloaddition of "66" bonds are car- 
ried out. To inspect the range of the exchange interac- 
tion, similar calculations are performed for the pair of 
CgoH radical adducts occupying the opposite (next near- 
est) vertexes of the tetragonal tetramer (C6o)4H2, which 
is a doped fragment of the 2D polymerized tetragonal 
phase of Ceo*^ 

The calculations have been carried out in the frame- 
work of the density functional theory (DFT) and ab 
initio Hartree-Fock (HF) methods. In the DFT calcu- 
lations, the hybrid functional of Becke (BSLYP)??!?^!?^ 
was employed, which includes the gradient-corrected ex- 
change and correlation functionals along with the exact 
exchange. The HF calculations were carried out using 
the PC GAMESS version^ of the GAMESS (US) QC 
packages^ The DFT calculations were carried out us- 
ing the PC GAMESSS and the Gaussian 03 suite of 
programs,-95. We exploited the spin-unrestricted method 
for both singlet and triplet spin states. The Gaussian 
basis sets employed are 3-21G and 6-31G*. The energy 
gradient convergence tolerance was less than 10^** A.U. 

Recently, the B3LYP method has been successively ap- 
plied to solids i^SiSLSSiSMSdi A significant improvement 
over the LDA results for electronic, structural and vibra- 
tional properties for some semiconductors and insulators, 
has been achieved by this method In addition, the 
B3LYP improves the magnetic moments and energy gaps 
and correctly predicts the ground state for some antifer- 
romagnetic insulators. ^^■^'^■^^ For the random impurity 
distribution case, the cluster approach is more appropri- 
ate then methods developed for periodic systems. 

The paper is organized as follows. In Section ^ we 
consider a dimer as the smallest nontrivial fragment of 
the lattice. We study the charge and spin distributions 
in the doubly charged dimer (Ceo)"*^^ and calculate the ex- 
change interaction of the spins localized on the fullerenes 
in the dimer. In Section IlIII we consider spin properties 
of dimers "doped" with radicals. To estimate the range 
of the exchange, we consider in Section IlIII a cluster with 
four fullerenes and calculate the next nearest neighbor 
exchange interaction. In the last section, we discuss the 
results and estimate the Curie temperature in our model. 
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II. CHARGED DIMER: THE EXCHANGE 
INTERACTION 

The rhombohedral (Rh) and tetragonal (Tg) polymer- 
ized Ceo crystals are built of Cgo layers in which fuUerene 
molecules are connected by [2+2] cycloaddition of "66" 
bonds. ^^'^^ The smallest lattice fragment is two adjacent 
Cgo molecules, a dimer, and we begin our study with the 
analysis of this simplest system in its neutral and charged 
states. 

It is known from the literature that the ground state 
of a free fuUerene dimer radical adduct i?-C6o-C6o-^ is 
a singly bonded isomer with the radicals R placed in spe- 
cific positions. ■^'^ The transition from the [2-1-2] double 
bond isomer to the singly bonded one requires break- 
ing of one of the bonds and a rotation of the buckyball 
around the remaining bond. This process may occur in 
a metastable dimer phase^S, but in a polymeric fuUcrcne 
network, where the molecules are tightly bound with the 
nearest neighbors, the rotation costs the bending and 
torsion energy and is, therefore, unfeasible. Thus, we 
assume that the charge and spin transfer from impuri- 
ties to fuUerene molecules take place without changing 
the bonding type. Therefore, we study a doubly charged 
[2-1-2] cycloadduct isomer (€60)2^ as a fragment of the 
doped polymerized fuUerene lattice. 

To obtain a quantitative information on the spin 
and bonding configuration, we have performed spin- 
unrestricted B3LYP (UB3LYP) calculations of the elec- 
tronic structure, the total energy, and optimized geome- 
try for the charged (€60)2^ dimer in the state with the 
total spin 5 = (singlet) and S — 1 (triplet). 

To check our methods, we began our calculation with 
the case of a neutral dimer. In agreement with earlier 
results?^ we found the [2-|-2]-cycloadduct isomer singlet 
with D2h symmetry to be the lowest energy state of 
(C6o)2- A triplet state of the doubly bonded neutral 
dimer lies 2.1 eV higher in energy. In the singlet state, 
the doubly charged [2-1-2] cycloadduct isomer has D2h 
symmetrysM The symmetry of the triplet states has not 
been studied before. We have considered the dimer in 
triplet state possessing one of the two simplest symmetry 
elements: inversion, Ci, symmetry and the mirror, Cg, 
symmetry with the plane of symmetry transverse to the 
dimer axis. The B3LYP calculations give identical results 
for the triplet states of both Ci and Cg isomers, and show 
that the ground state of [2-1-2] cycloadduct (050)2^ dimer 
is the D2h isomer in the triplet spin state. 

Listed in Tables U and ^ the results of the 
B3LYP/6-31G* calculations for both negatively charged 
(C6o)2~ ^-iid positively (050)2''' charged dimers show that 
the ground state of a doubly charged [2-t-2]-cycloadduct 
isomer is the triplet spin state of the D2h symmetry. 
It is interesting that the properties of the negatively 
( "electron doped" ) and positively ( "hole doped" ) charged 
dimers are rather close, as seen from the comparison of 
Tables III and ITU Also, one sees from Tables III Ull that the 
results obtained at 3-21G and 6-31G* levels are in close 



TABLE L The total energies, lengths and orders of the in- 
terfuUerene bonds for T)2h isomer (C6o)2~ in the singlet and 
triplet states for the optimized structure calculated by the 
use of the UB3LYP hybrid DFT method with the 3-21G and 
6-31G* basis sets employed 



basis set 


3-21G 6-31G* 


spin state 


singlet triplet singlet triplet 


total energy (A.U.) 
bond length (A) 
bond order 


-4547.166 -4547.175 -4572.430 -4572.442 
1.593 1.597 1.597 1.593 
0.799 0.796 0.859 0.864 


TABLE IL The total energies, lengths, and orders of the 
interfuUerene bonds for T)2h isomer (050)2''' singlet and 
triplet states for the optimized structure calculated by the use 
of the UB3LYP hybrid DFT method with the 3-21G and 6- 
31G* basis sets employed. 


basis set 


3-21G 6-31G* 


spin state 


singlet triplet singlet triplet 


total energy (A.U.) 
bond length (A) 
bond order 


-4546.434 -4546.442 -4571.765 -4571.772 
1.584 1.602 1.581 1.597 
0.806 0.791 0.874 0.859 



agreement. 

To identify the nature of the stationary points on the 
potential energy surface (a true minimum or a saddle 
point), we calculated the vibrational frequencies for the 
neutral singlet -D2/1 (C6o)2 isomer and for both the sin- 
glet and triplet -D2/1 (050)2^ isomers. The absence of 
imaginary vibrational frequencies indicates that all the 
stationary points are minima. 

A robust feature seen from the ab initio dimer calcula- 
tions is that the bridging bonds of the [2-1-2] cycloadduct 
are not affected by doping and are almost impenetra- 
ble for the spin. In other words, each fuUerene in a 
charged dimer possess a well defined spin localized on the 
fuUerenes. Hence, one can conveniently describe the in- 
teracting pair in terms of the Heisenberg Hamiltonian,— 

H = -2JSi ■ S2 , (2) 

where Si (5'^ = |) is the spin vector on the i-th site of 
the dimer, and J is the corresponding exchange integral. 
The latter can be presented as 

J={En-E^{)/2, (3) 

where Ei^i and are the energies of the dimer in the 
singlet and triplet states, respectively. 

We calculate the exchange integral J, from the singlet- 
triplet splitting Eq. O having taken the corresponding 
energies from Tables and The results are listed in 
Table Hill For all the calculation methods we obtained a 
positive exchange integral, i.e., the interaction between 
the spins localized on the adjacent fuUerenes is ferromag- 
netic. Both spin-unrestricted and spin-restricted B3LYP 
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TABLE III: Effective exchange integral J = {E-]i - Eff) for 
pair of adjacent ions in polymerized fuUerene. The ener- 
gies El I and Eii are total energies of low- lying singlet and 
triplet states for corresponding [2+2] cycloadduct charged 
D2h isomers 





J (eV) 


method 


(C6o)2 


(C6o)r 


AMI 


0.52 


0.45 


B3LYP/3-21G 


0.12 


0.10 


B3LYP/6-31G* 


0.16 " 


0.10 



"The ab initio Hartree-Fock/6-31G* calculations give in this case 
the value equal to 0.72 eV 



calculations give the same value of the singlet state en- 
ergy and, therefore, the singlet spin state is not "antifer- 
romagnetic" by its nature. Note that both the semiem- 
pirical AMI and ab initio Hartree-Fock methods, where 
correlation effects are ignored, noticeably overestimate J. 
The B3LYP hybrid functional calculations, which take 
the electron correlations into account along with the ex- 
act exchange, are more accurate and reliable; they give 
far lower values of the exchange integral: J=0.16 eV for 
(C6o)2" and 0.10 eV for (€60)2+. 

This is the key result supporting our model of fer- 
romagnetic fuUerenes: in all the studied cases of the 
[2-1-2] cycloadducts (€60)2^ , the exchange integral is in- 
variably positive and the exchange interaction is rather 
strong. We emphasize that the ferromagnetic exchange 
is an intrinsic property of polymerized [2-1-2] cycloadduct 
fuUerenes. 

We note also that the spin intercage interaction is more 
sensitive to the relative orientation of the fuUerenes then 
to the intercage separation: We see from ab initio calcu- 
lations that for the relative orientation which corresponds 
to a singly bonded dimer, the separation is almost the 
same as in the [2-1-2] cycloadduct case, but the ground 
state for the singly bonded isomer (C6o)2^^ is different: It 
is spin singlet separated from the triplet state by the gap 
~ 0.7 eV at the B3LYP/3-21G level. 



III. EXCHANGE INTERACTION FOR 
RADICAL ADDUCTS Ceoi? 

Most important question is the availability of impu- 
rities with the wanted property to create spins local- 
ized on fuUerenes in the polymerized matrix. Two types 
of doping can be envisaged, depending on the impurity 
character. The first one is a charged complex formed 
in the fuUerene matrix, comprising a fully ionized impu- 
rity bound primarily by the Madelung electrostatic forces 
(ionic limit). This is the case, for example, in TDAE- 
Cgo, as well as for the alkali metal doping. Second, an 
impurity atom or molecule forms Cgo radical adduct cova- 
lently bound with a carbon atom of the Cgo cage (covalent 
limit). This is the case, for example, for hydrogen, flu- 
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FIG. 1: The charge (a) and spin (b) density spatial dis- 
tribution (MuUiken atomic charges and spins) of the doped 
fuUerenes: CeoH, CeoF. Black (grey) corresponds to positive 
(negative) charge or up (down) spins. The volume of the 
spheres is proportional to the absolute value of the corre- 
sponding variable on a given atom. 



orine, hydroxyl group OH, amino group NH2, or methyl 
group CH3. 

We will focus our attention on the covalently bounded 
radicals (ligands) and especially on hydrogen. Particu- 
lar interest in hydrogen doping stems from the fact that 
hydrogen, a donor, is always present in fuUerene solids 
in noticeable amounts. Hydrogen was detected in rather 
large concentration (about one hydrogen atom per six 
fuUerene molecules) in ferromagnetic samples of pressure- 
polymerized fullereneii24 Fluorine is of general interest for 
it is one of the very few atoms that has stronger acceptor 
properties than Cgo. Hydroxyl, amino and methyl groups 
exemplify simplest molecular radicals. 

We begin with the presentation of our results concern- 
ing the electronic structure the above radical adducts. 
The spatial distribution of the MuUiken atomic charge 
and atomic spin for CeoH and CgoF calculated at 
B3LYP/3-21G level is presented in Fig. □ Atomic Mul- 
liken charges (the sum of the atomic MuUiken charges in 
the case of fuUerene molecule), reflect the nature of the 
chemical bond, i.e., the degree of hybridization of the va- 
lence orbitals of the fuUerene and impurity atoms, rather 
than the values of atomic net charges. The net charges 
are closely related to the spin density distribution. For 
CeoH and CgoF, the B3LYP1/3-21G calculations show 
that the MuUiken charges of H and F are equal to 0.27 
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a b 

FIG. 2: The charge (a) and spin (b) density spatial dis- 
tribution (MuUiken atomic charges and spins) of the doped 
CA dimers (C6oH)2, (C60F)2, (C6o)2HF. Black (gray) corre- 
sponds to positive (negative) charge or up (down) spins. The 
volume of the spheres is proportional to the absolute value of 
the corresponding variable on a given atom. 



and -0.25, correspondingly, whereas the atomic spin pop- 
ulations for both H and F atoms are equal to 0.042. The 
latter values imply that the net charges for H and F are 
equal to ±0.958, respectively, and corresponding Ceo net 
charges are equal to =f0.958. For comparison, sums of 
atomic spin populations for hydroxyl group OH, amino 
group NH2, and methyl group CH3 are equal to 0.043, 
0.072, and 0.042 respectively, that is much less than unity 
as for H and F, whereas corresponding sums of MuUiken 
atomic charges are equal to -0.16, -0.03 and 0.15. 

To estimate the exchange interaction between the 
spins on neighboring fuUerenes in the matrix, we have 
performed B3LYP/3-21G calculations for the cases of 
doped [2-1-2] cycloadduct dimers H-C6o=C6o-H and 
F-C6o=C6o-F, OH-C6o=C6o-OH, NH2-C6o=C6o-NH2, 
CH3-C6o=C6o-CH3, as well as for a combination of 
donor and acceptor H-C6o=C6o-F and H-Cgo^Cgo-OH. 
Several configurations differing in the initial position of 
the ligands relative to the fuUerene molecules have been 
considered. For all of the configurations, the spins on the 
neighboring fullerenes are parallel in the ground state, 
forming a triplet, and the energy differences E^i — E-^^ 
are of the order of several tenths of eV. The triplet-singlet 
splittings for the configurations shown in Fig.|21are listed 
in Table ITVI 

The spatial distributions of the charge and spin density 
for doped dimers i?-C6o=C6o--R' (with R, R' — H, F) are 
shown in Fig. |21 The spin density. Fig. [2Ib), is spread 
across the buckyballs, repeating the net charge density 
profile. Comparing Fig. [3 and Fig. one concludes that 
in the dimer case the spin distribution is nearly the same 



TABLE IV: Total energies for the singlet. En, and triplet, 
E^-, states, and energy gain of spin polarized state, E^ — E^-, 
calculated at the UB3LYP/3-21G level for a number of doped 
[2-1-2] cycloadduct 7i— Ceo =C6o — R' dimers and for planar 
tetragonal tetramer (C6o)4ll2. 



doped isomer 


E-\i 


Ell 






(A.U.) 


(A.U.) 


(eV) 


H~C6o=C6o— H 


-4548.182 


-4548.203 


0.58 " 


F-C6o=C6o-F 


-4745.555 


-4745.578 


0.61 


H— C6o=C6o— F 


-4646.869 


-4646.890 


0.56 


OH— Ceo^Ceo— OH 


-4697.751 


-4697.773 


0.59 


H-C6o=C6o-OH 


-4622.967 


-4622.988 


0.58 


NH2-C60=C60-NH2 


-4658.246 


-4658.268 


0.58 


CH3-C60=C60-CH3 


-4626.393 


4626.414 


0.57 


(C6o)4H2 


-9095.234 


-9095.259 


0.68 



"The B3LYP/6-31G* calculations give in this case the value equal 
to 0.57 eV 



as for an isolated Cgo^. Note that the spin density dis- 
tribution is insensitive to the sign of the fuUerene charge. 
Also, the spin and net charge distributions are insensitive 
to the choice of the radical: Our numerics shows that for 
i?-C6o=C6o-i?, R= OH, NH2, CH3, the spatial distribu- 
tions of the net charge and spin density over the fuUerene 
molecules are very similar to those in Fig. |21 

Remarkably, the spin density is zero on the bonds 
bridging the molecules, so that one can assign a spin to 
each of the buckyballs. This corresponds to the picture 
where the spin of an individual charged (doped) molecule 
is well defined, and Ceo molecules play the role of the 
sites on which spins reside. From Table ITVI the nearest 
neighbor exchange interaction Eq. |(2Jl is rather strong, 
J « 0.3 eV. 

To estimate the spatial range of the exchange inter- 
action, we consider a pair of fuUerene hydrogen adducts 
CgoH placed in the next nearest neighbors positions. For 
this, we compute properties of planar tetragonal tetramer 
(C6o)4H2 (Fig.O with doped Cgo at the opposite corners 
of the tetramer. The UB3LYP/3-21G calculations of the 
electronic structure and the total energy have been car- 
ried out for optimized geometry of the singlet and triplet 
spin states of the cluster. The singlet-triplet splittings 
are presented in Table IIVI The charge and spin den- 
sity spatial distributions for the tetramer (C6o)4H2 dis- 
played in Fig. O are very similar to that for the dimer 
in Fig. 121 The singlet-triplet splitting is positive (ferro- 
magnetic) and close to that in a dimer. Therefore, we 
observe that the exchange interaction between fuUerene 
radical adducts in a 2D layer does not fall within the 
two first coordination spheres. Moreover, the indirect ex- 
change interaction mediated by fuUerene matrix is some- 
what stronger than direct one. Thus, it is very likely 
that the spin interaction extends far beyond the nearest 
neighbors, keeping its sign and magnitude. Of course, 
additional calculations of spin interaction between dis- 
tant fuUerene radical adducts in a large fragments of the 
lattice are needed to support this conjecture. 
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FIG. 3: The charge (a) and spin (b) distribution in the triplet 
spin state of the doped CA tetramer (C6o)4H2. Black (gray) 
corresponds to positive (negative) charge or up (down) spins. 
The volume of the spheres is proportional to the absolute value 
of the corresponding variable on a given atom. 

In conclusion of this section, we emphasize that the fer- 
romagnetic exchange is an intrinsic property of polymer- 
ized [2-1-2] cycloadduct fullerenes rather than the radical 
impurity that makes the fuUerene magnetically active. 



IV. DISCUSSION AND CONCLUSIONS 

Based on the presented ab initio calculations, the fol- 
lowing physical picture emerges. Chemically pure poly- 
merized fuUerene lattice, rhombohedral or tetragonal , is 
a diamagnet with zero spin on the sites of the Cgo ma- 
trix. Doping with radicals initiates reaction Eq. and 
some of the cites of the matrix become paramagnetic as 
the result of the spin (and charge) transfer from the rad- 
ical attached to the corresponding fuUerene. The calcu- 
lations of the spin distribution on the doped dimers and 
tetramers show that the spin transferred from the radical 
is smeared over the buckyballs surface, but it avoids the 
[2-1-2] four-membered rings connecting Cgo's, so that the 
transferred spin is well localized on the fullerenes. One 
comes to the picture of immobile spins occupying sites of 
the fuUerene matrix. 

The interaction of the localized spins can be described 
by the standard Heisenberg Hamiltonian. In accordance 
with our calculations, the Heisenberg exchange integral 
J is ferromagnetic and rather strong, J — 0.1 ^ 0.3 eV. 
The cluster calculations in Section Hill give the evidence 
that the ferromagnetic exchange extends at least to the 
next nearest neighbor and, probably further. Polymer- 
ized fullerenes, rhombohedral and tetragonal, are quasi 
two-dimensional crystals, and it should be noted that 
the strong exchange found in our calculations refers only 
to the in-plane interaction. We assume the interplane 
exchange interaction J' to be ferromagnetic and weak, 
similar to that in the low Tc Ceo-TDAE-compounds. Its 
exact value is not of primary importance, as discussed 
below. 



In our model, a ferromagnetic transition occurs in a 
solid solution CeoRx where the impurity (dopants) con- 
centration X is small. For small x, the spatial distribution 
of the dopant is random so that the sites of the fuUerene 
matrix become spin active in a random fashion. Con- 
sequently, doped polymerized fullerenes are expected to 
belong to the class of disordered ferromagnets. Proper- 
ties of a disordered magnet are commonly discussed in 
the framework of percolation theory (see Ref. and ref- 
erences therein). A detailed analysis is beyond the scope 
of the paper, and we limit ourselves to few remarks. 

At zero temperature, any two spins are aligned if they 
are within the range of the exchange interaction. The 
long-range magnetic order, when the aligned spins belong 
to the infinite cluster, is established provided the spin 
concentration, that is the impurity radical concentration, 
exceeds the percolation threshold. For a two-dimensional 
square lattice with only the nearest neighbors interaction, 
the critical concentration ~ 0.59. If the interac- 

tion extends beyond the nearest neighbors, the critical 
concentration scales as Ti~^ with the radius of the ex- 
change interaction TZ (measured in units of the lattice 
constant). Our calculations show that the interaction 
is not limited to the nearest neighbors, and TZ is not less 
than 2V2 and then Xc < 0.07. In accordance with this es- 
timate, a single impurity radical for 14, or perhaps more, 
Ceo molecules sufUces to transform a diamagnetic poly- 
merized fuUerene into a ferromagnet. 

The assumption of the percolation theory about a ran- 
dom distribution of impurities may be invalid for certain 
synthesis conditions, when segregation of defects may oc- 
cur. In this case, the formation of ferromagnetic islands 
is expected, with or without bulk magnetic order. 

Evaluation of the Curie temperature, Tc, is hindered 
by the presence of disorder and a quasi two-dimensional 
nature of the polymerized fullerenes. In accordance with 
the Mermin- Wagner theorem^, two-dimensional (2D) 
isotropic magnets exhibit a long-range order only at 
r = K for any value of the in-plane exchange inte- 
gral J. In quasi-two dimensional layered compounds, the 
magnetic order establishes at a finite temperature Tc due 
to the inter-plane exchange J'. It is well-known^* that in 
the limit of small J', the Curie temperature can be esti- 
mated as Tc ~ J/ In -p- ; the estimate should be valid for 
our case of a disordered system unless very close to the 
percolation threshold. In accordance with our ab initio 
calculations, the in-plane exchange interaction is in the 
range J = 0.1 0.3eV = 1200 3500 K . Seeing that the 
dependence of Tq on J' is only logarithmic, the Curie 
temperature Tq is several times (but not several orders) 
less than J, i.e., in the range 300 -f- lOOOi^. 

These estimates, although crude, show that the model 
has the potential to give the interpretation to the 
experiments^i^*S4L2424i& where magnetism of polymerized 
fullerenes was observed at room temperatures and above. 
These calculation show also that a low magnetization, 
which is controlled in our model by the concentration of 
the radicals, is compatible with a high Curie temperature. 
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The presence of any radicals in the studied samples 
of polymerized fuUerenes is an open question. The only 
exception is hydrogen, which was detected in a rather 
large amount (one hydrogen atom per six fuUerenes) in 
some polymerized fuUerenes exhibiting high-temperature 
ferromagnetismi^. Some caution is needed at this point. 
Hydrogen is known to have the tendency towards the for- 
mation of the radical adducts Ceo Hn with even number 
n. Our calculations show that the ground state of Ceo-ffn 
complex is spin singlet for n = 2, and most likely this 
conclusion holds for higher even n's. For this reason, the 
presence of hydrogen in large amounts does not guarantee 
magnetism. However, if we are concerned with a dilute 
solid solution of hydrogen, where rare hydrogen atoms 
are dispersed chaotically in the fuUerene matrix, the for- 
mation of CeoH complexes becomes feasible during the 
high-temperature synthesis. Formally a metastable state, 
the system of the CgoH complexes quenched upon cool- 
ing, is practically stable, being separated from the true 
energy minimum by a high energy barrier (of the order of 
2 eV per pair of hydrogen atoms in accordance with our 
numerical calculations). In our scenario, the quenched 
network of the magnetically active CeoH sites becomes 
ferromagnetic below Tc- The formation of the network 
is obviously sensitive to details of the synthesis condi- 
tions, and from this point of view, a poor reproducibility 
of the magnetic fuUerene synthesis procedure can be un- 
derstood. 

Besides hydrogen, the presence of fluorine, hydroxyl 
group OH, amino group NH2, and methyl group CH3, 
is favorable, in accordance with our ab initio calculation. 



for the high-temperature ferromagnetic phase of polymer- 
ized fuUerenes. Since these materials have not yet been 
synthesized and measured, this is a theoretical possibility 
waiting for experimental verification. 

As has been already mentioned, a similar mechanism 
of fcrromagnctism is responsible for a low temperature 
transition of the Cgo fuUerene intercalated with the donor 
TDAE molecules. The exchange interaction is much 
weaker in this case^ because the mutual orientation of 
the adjacent buckyballs differs from that in the polymeric 
phases. 

In conclusion, we have considered a mechanism of fcr- 
romagnctism in polymerized fuUerenes, where diamag- 
netic Geo molecules transform into stable paramagnetic 
species, ions C^q or fuUerene radical adducts CqqR, and 
become magnetically active due to the charge and spin 
transfer from radical impurities. The model is supported 
by ab initio calculations, the main result of which is that 
in the [2-|-2]-cycloadduct polymerized phase the effective 
exchange interaction between the paramagnetic species 
is ferromagnetic and strong enough to account for the 
high temperature ferromagnetism observed in recent ex- 
periments on polymerized fuUerenes. The model predicts 
ferromagnetism with the high Curie temperature in the 
polymerized fuUerenes doped with the radicals like hy- 
drogen, fluorine, hydroxyl group, amino group, or methyl 
group. 
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